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ABSTRACT 

Transformation of disks into spheroids via mergers is a well-accepted element of galaxy 
formation models. However, recent simulations have shown that bulge formation is suppressed 
in increasingly gas-rich mergers. We investigate the global implications of these results in a 
cosmological framework, using independent approaches: empirical halo-occupation models 
(where galaxies are populated in halos according to observations) and semi-analytic models. 
In both, ignoring the effects of gas in mergers leads to the over-production of spheroids: 
low and intermediate-mass galaxies are predicted to be bulge-dominated (B/T ~ 0.5 at 
< 10 10 Mq, with almost no 'bulgeless' systems), even if they have avoided major mergers. 
Including the different physical behavior of gas in mergers immediately leads to a dramatic 
change: bulge formation is suppressed in low-mass galaxies, observed to be gas-rich (giving 
B/T ~ 0.1 at < 1 10 Mq, with a number of bulgeless galaxies in good agreement with obser- 
vations). Simulations and analytic models which neglect the similarity-breaking behavior of 
gas have difficulty reproducing the strong observed morphology-mass relation. However, the 
observed dependence of gas fractions on mass, combined with suppression of bulge formation 
in gas-rich mergers, naturally leads to the observed trends. Discrepancies between observa- 
tions and models that ignore the role of gas increase with redshift; in models that treat gas 
properly, galaxies are predicted to be less bulge-dominated at high redshifts, in agreement 
with the observations. We discuss implications for the global bulge mass density and future 
observational tests. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — galaxies: spiral 
— cosmology: theory 



1 INTRODUCTION 

In a ACDM cosm ology, structure grows hierarchically (e.g. 
IWhite & ReesI [l978h . making mergers and interactions between 
galaxies an essential and inescapable process in galaxy forma- 
tion. Indeed, mergers are widel y believed to be responsible for 
the morphologies of spheroids ( lToomrelll977l) . and observations 



find recent merger remnants in considerable abundance in the lo 
cal universe (Schweizer 1982; Lake & Dressier 1986; Dovon et al. 
1994l;IShier& Fischer 1998; I James et al. 1999; Genzel et alfeOOll; 



■i; 

6) 



Tacconi et alj|2002t iRothberg & JosephH2004l ; iDasvra et allHooi 
as well as e.g. faint shells and tidal features common around ap- 
paren tly "normal" galaxies ( Mali n & CarteJl98dl : ISchweizeJl98d. 
Il996h. which are thought to be signatures of galaxy collisions (e.g. 
Hernq uist & Ouinnlll988l : lHernquist & Spergelll 1992b . 
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This has led to the concern that there may be too many 
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mergers to explain the survival and abundance of galactic disks 
in the c ontext of our present understanding of g alaxy formation 
(see e.g. IStewart et alj|2008bl : IPurcell et al] 120091 and references 
therein). Indeed, the difficulty of forming very late type disks 
with low B/D is a well-known problem in cosmological simula- 
tions, which tend to consistently overpredict the bulge fractions 
of low-mass galaxies. There has been considerable effort to ad- 
dress the cause of these difficulties, with a great deal of focus on 
improving resolution (in order to properly resolve the bulge and 
disk components) and introducing feedback prescriptions to pre- 
vent artificial disk over-coo ling, clumping, fragmentation, and an- 
gular momentum transport dWeil et alll998l;ISommer-Larsen et al.l 
19991 l2003l:lThacker & Couchmanll2000l.l200ll:lAbadi et alj|2003l; 



Governato et alj|2004 2007 . Robertson et al.ll2004IOkamoto et al.l 



20051 ; IScannapieco et alj|2008h . Similar problems are encountered 
in semi-empirical models, that use observational constraints to pop- 
ulate a halo-occupation fra mework, and evol ve this forward in time 
to follow mergers (see e.g. lStewart et al.f2009l) . A number of semi- 
analytic models, which predict galaxy properties from simple phys- 
ical recipes and adopt simplified prescriptions for behavior in merg- 
ers, allowing for the efficiency of merger-induced bulge formation 
to be adjuste d, still find the same difficulty even with these degrees 
of freedom dSomerville et alj 1200 ll, l2008ij; iKhochfar & Burkertl 
l2003l : lGranato et all2004ICroton et alj|2006h . 

More importantly, the problem can be seen empirically with- 
out reference to these models. Observations find that ~ 5 — 10% 
of low or intermediate-mass galaxies (M* < 10 10 Mq at redshifts 
z ~ 0.2 — 1.2) are in mergers - strongly morphologically disturbed 
or in "major" similar-ma ss pairs about to merge (at small scales and 
small relative velocitie s) jBridge et al. 2007, 2009; Kartaltepe et al 



20071: IConselice et alj|2009l. 120081 ; IJogee et alj [20091 ; lLotz et al 
2008b; iLin et al.l l2008h . This is equal to or higher than the pre 



dieted rates from the models above - the problem does not seem to 
be that ACDM predicts "too many" mergers. The problem is that, 
even without correcting for the short expected duty cycles/lifetimes 
of this merger phase, these fractions are already higher than the 
observed fraction of bulge-dominated galaxies at these masses 
dAllen et al.l2006l;lBalcells et al.l2007bl ; lDomi'nguez-Palmero et alj 
120081 ; IWeinzirl et al.ll2009l) . Invoking subsequent accretion to "re- 
build" disks will not eliminate massive bulges formed by these 
events, nor can it operate in a shorter time than the duration 
of the merger itself, and correcting for the expected duty cy- 
cles/lifetimes of these events, a large fraction of these low-mass 
galaxies s hould have experienced a violent event in the last 
~ lO Gyr dHopkins et a"Hl2007l ; lLotz et alj|2008 j ; IConselice etail 
2009). Moreover, although there may be some mass-dependence in 
observed merger fractions, it i s weak: a factor ~ 2 increase with 
mass from ~ 10 9 - 10" M p dZheng et alj|2007 l: IConselic e et all 
2008MPatton & A tfield 2008 ! :lBundv et alj2009l:lBridge et alj2009l ; 



Lopez-Saniuan et al. 20091 : Domingue et al. 20091) . Over the same 
range, there is an order-of-magnitude or more change in the fraction 
of bulge-dominated galaxies. As a result, uniformly increasing or 
decreasing the "efficiency" of merger-induced bulg e formation is 
insufficient to explain the observed trends (see e.g. IStewart et al.l 
2009). In short, the observed mergers at low and intermediate 
masses must somehow be inefficient at making bulge s, whereas 
those at high masses must be efficient ( Koda et al. 2009). However, 
dissipationless systems with high or low stellar and/or virial masses 
behave exactly the same in actual A-body simulations of mergers. 
They must do so — if gas is not allowed to change the dynam- 
ics of the merger, then the only force is gravity and the system is 
scale-free (a merger of two 10 9 Mq systems is exactly equivalent to 



the merger of two 10 1_ Mq systems). Because gravitational inter- 
actions of dissipationless systems are scale-free, it is not possible 
to design a physical prescription that does not treat gas differently 
from stars (or where the merger efficiency is just a function of mass 
ratio and orbital parameters) that can simultaneously match the 
observed merger fractions and B/T distributions at low and high 
galaxy masses. Some other controlling parameter or physics must 

wor k to suppress bulge formation in low- mass galaxies. 

ISpringel & Hernquistl |2005) and iRobertson et al.1 d2006al) 
suggested that this parameter may be galaxy gas fraction^ they 
showed, in idealized merger simulations, that even major mergers 
of sufficiently gas-dominated disks could produce remnants with 
large di sks. This has since been bor ne out in cosmological simu- 
lations dGovernato et alj|200ll2008l) . and may be seen in observa- 
tions .jrfJiigJw^cMvifufe merg- 
ers dGenzeletalJl2008l: [ Shapiro et all 120081: iKassin et alj 120071: 
van Dokkum et alj|2004 IPuech et alj|2008l : Ivan S tarkenbur g et all 
20081 : IRobertson & BullocU l2008a). The growing consensus that 



more realistic treatments of star formation and supernovae feed- 
back - which suppress efficient early star formation - play a crit- 
ical role in enabling disk formation in simulations (see references 
above) is connected with this insight into the link between gas frac- 
tion and disk survival. Because low-mass disks are observed to be 
more gas-rich, this may potentially provide the desired breaking of 
self-sim ilarity in dissip a tionless merger histories. 

In iHopkins et all d2009bh (hereafter H09), we developed 
a physical model to explain this be havior, and used a large 
suite of hydrodynamic simulations (see lCox et al . 2006b ai. 120081 : 
IRobertson et al.l2006d lbT) to show that the suppression of bulge for- 
mation in mergers of more gas-rich disks can be predicted analyt- 
ically, and is dynamically inevitable owing to fundamental prop- 
erties of gravitational physics and collisional systems. Stars, be- 
ing collisionless, violently relax in mergers, as they can mix and 
randomize their orbits by scattering off the fluctuating potential, 
even without net torques. Gas, being collisional, cannot violently 
relax, but must have its angular momentum torqued away in order 
to dissipate and build a bulge by forming stars in a central star- 
burst. In a merger, this torque is primarily internal, from stars in 
the same disk: the passages and merger of the secondary induce 
a non-axisymmetric distortion in the primary stellar and gas disk. 
The stellar distortion (in a trailing resonance and close proximity 
to the gas response) efficiently removes angular momentum from 
the gas and allows it to fal l to the center and form a bulge (see also 
iBarnes & Hernquistl 19961) . Other sources of torque (hydrodynamic 
shocks or the gravity of the secondary galaxy) represent less effi- 
cient, higher-order effects (and often tend to spin up the disk). 

In the limit, then, of extremely gas-rich disks, there is no sig- 
nificant stellar feature available to torque the cold gas in the disk; 
gas will rapidly re-form a cold disk in the background of the re- 
laxing stellar potential. Because this is a purely gravitational pro- 
cess, the physics and consequences for fixed initial conditions are 
independent of uncertain physics associated with "feedback" from 
either star formation and supernovae or accreting black holes (see 
H09). Feedback is important for setting these initial conditions - 
over cosmological timescales, it can suppress star formation ef- 
ficiencies and lead to retention of large gas reservoirs, and move 



1 In this paper, we use the term "gas fraction" to refer specifically to the 
cooled gas fraction in the baryonic disk, i.e. (mm + «iH 2 )/( m Hl + "JH 2 + 
m» ), where i»hi, oih 2 > and m * me the atomic, molecular, and stellar masses, 
respectively. 
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that gas to larger radii; but for a given set of gas properties at the 
time of the merger, H09 outline a few simple equations motivated 
by well-understood gravitational physics that can robustly estimate 
how much of the gas will preserve its angular momentum content 
in a given merger. 

We endeavor here to study the consequences of these detailed 
studies of disk survival as a function of merger gas content, for the 
global demographics of bulges and disks as a function of galaxy 
mass and redshift. To do so, we utilize both semi-analytic and 
semi-empirical cosmological models. Semi-analytic models predict 
galaxy properties from first principles via physical recipes; while 
the latest generation of models reproduce many observational quan- 
tities quite well, there remain significant discrepancies and uncer- 
tainties. We therefore compare to a semi-empirical model, where 
we assign galaxy properties to dark matter (sub)-halos in order to 
match observational constraints. We combine these with the de- 
tailed predictions from simulations presented in H09, to make pre- 
dictions for bulge and disk populations. We contrast these predic- 
tions with the results from implementing the "traditional" assump- 
tions regarding bulge formation, which ignore the unique role of 
gas as outlined above. 

In §[2] we describe our empirical approach for tracking stellar 
and gas masses within the framework of our cosmological merger 
trees, and briefly summarize our semi-analytic approach. In § [3] 
we outline the key predictions of our gas fraction-dependent bulge- 
formation recipe and the differences from those of the recipes gen- 
erally adopted in the past. We summarize and discuss other impli- 
cations in §[4] 

Throughout, we assume a flat cosmology with matter and dark 
energy densities Q.m = 0.3 and Q.a = 0.7, respectively, and Hub- 
ble constant Ho = 70kms -1 Mpc~', but var ying these parameters 
within the observationally constrained limits dKomatsu et al.l2009l) 
makes little difference. 



2 THE MODEL 

In order to track merger histories with as few assumptions as pos- 
sible, we construct the following semi-empirical model, motivated 
by the halo occupation framework. Essentially, we assume galax- 
ies obey observational constraints on disk masses and gas fractions, 
and then predict the properties of merger remnants. 

We follow a Monte Carlo population of hal os from cosmo- 
logica l simulations, with merger/growth trees from Fakhouri &Ma 
(2008). We populate every halo with a 

gaiaxJT 

on the 

empirically constra i ned s tellar mass-halo mass relation from 
IConrov & WechsleJ [2009). As the halo grows, the galaxy is as- 
signed a star formation rate and net gas accretion rate (inflow mi- 
nus outflow) corresponding to the track that would be obtained by 
strictly enforcing that galax ies at all times live on t he same stellar 
mass-halo mass relation dConrov & Wechslerj[2 009 ) and gas mass- 
stellar mass rel ation (fitted from ob servations as a function of mass 
and redshift in IStewart et al]|2009h . In other words, using the em- 
pirical fits to Mg a i(Mhaio) and / ga s(M ga i), we force galaxies to live 
on these relations at all times assuming gas is gained or lost, and 
stars are formed, exactly as needed to stay on the empirically fit- 
ted correlations, but without a prior on how much of each is in 



2 Starting at some minimum M^\ ~ 10 10 Mq, although this choice makes 
little difference. 



a bulge or disk component. This amounts to the same prescrip- 
tion as assuming galaxies l ive on e.g. the obse rved stellar mass- 
star formation rate relation (No eske et alj|2007l) and tuning accre- 
tion rates to reproduce the observed dependence of gas fractions 
on stellar mass and redshift teell & de Jondl2000l;lMcGaughll2005l ; 
ICalura et alj|2008l ; IShaplev et al.ll2005l ; lErb et alj[2006t) . Observa- 
tional constraints are present up to redshifts z ~ 3, more than suf- 
ficient for our low-redshift predictions. Despite its simplicity, this 
effectively guarantees a match to various halo occupation statistics 
including stellar mass functions and the fraction of active/passive 
galaxie s as a function of m ass jYangetalJ|200l ;l Wein mann et al.l 
l2006allbllGe"rke et alj2007h . 

Together, these simple assumptions are sufficient to define 
a "background" galaxy population. There are degeneracies in the 
model - however, we are not claiming that this is unique nor that 
it contains any physics thus far. For our purposes, the precise con- 
struction of the empirical model is not important - our results are 
unchanged so long as the same gas fraction distributions are repro- 
duced as a function of galaxy mass and redshift. 

We illustrate this by comparison with a semi-analytic galaxy 
formation model dSomerville et al.|[2008al . hereafter S08). A sum- 
mary of the salient model features is presented in Appendix[B] Un- 
like the semi-empirical model above, the semi-analytic model at- 
tempts to predict the relevant quantities (e.g. disk masses and gas 
content) from first principles, given physical prescriptions for cool- 
ing, star formation, and feedback. Specifically, the model presented 
begins with similar halo merger trees, allowing for subhalo evolu- 
tion and stripping, and distinguishes between accretion in the "cold 
mode" in sub-L„ galaxy halos (where gas accreted into the halo 
reaches a galaxy on a dynamical time) and "hot mode" in massive 
halos (where the cooling occurs from a hot atmosphere, suppressed 
by AGN feedback). The model adopts typical recipes for star for- 
mation and supernova feedback, motivated by observations. 

The relevant physics are entirely contained in the treatment 
of mergers. When a merge£] occurs, we adopt a model for how 
much of the gas is drained of angular momentum and participates 
in a nuclear starburst, and how much of the stellar disk is violently 
relaxed — both components are removed from the disk component 
and added to the mass of the galaxy bulge. If orbital parameters — 
namely the relative inclination angles of the merging disks — are 
required, we simply draw them at random assuming an isotropic 
distribution of inclinations (allowing some moderate inclination 
bias makes no difference). We consider three models: 

(1) Gas Fraction-Dependent ("Full") Model: The prescrip- 
tions from H09, based on a detailed physical model for how 
disks are destroyed in mergers and calibrated to the results of 
hydrodynamic merger simulations. In a merger of mass ratio /i, a 
fraction ~ fi of the primary stellar disk is violently relaxed, and 
a fraction ~ (1 — / ga s) /J. of the gas loses its angular momentum. 
For the full details of these scalings as a function of mass ratio, 
orbital parameters, and gas fraction, we refer to H09, but they are 
briefly listed in the Appendix lAl the general sense is that angular 
momentum loss is suppressed by a factor ~ ( 1 — / gas ) . 

(2) Gas Fraction-Independent ("Dissipationless") Model: 

The same, but without accounting for the different behavior of 
gas and stars in the merger. We adopt the scaling from H09, 



3 We include all mergers above a minimum mass ratio fj, ~ 0.01, although 
our results are not sensitive to this limit, so long as it is small. 
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but calculate the angular momentum loss in the gas as if the 
disks were entirely stellar — equivalently, we ignore the fact that 
gas-richness changes the efficiency of bulge formation in mergers. 
In order for gas to behave in this manner it would have to be 
dissipationless and collisionless; but nevertheless this is a common 
(almost ubiquitous) assumption. The specific scalings are given in 
Appendix lAl 

(3) Simplified ("Threshold") Model: A "threshold" ("ma- 
jor" versus "minor") assumption, common in many semi-analytic 
models. In major mergers {p, > 0.3), the model assumes that 
the entire stellar disk is violently relaxed and the entire cold gas 
content participates in a burst. In minor mergers, the primary is 
unperturbed, the stars of the secondary are added to the bulge, and 
the secondary gas content is added to the disk. 



3 CONSEQUENCES 

3.1 Bulges and Disks at z = 

Figure Q] shows the predicted distribution of bulge-to-total (bulge 
stellar mass divided by total galaxy stellar mass) ratios as a func- 
tion of galaxy mass. Figure[2]co mpares these pr edictions to obser- 
vations of field spirals from Weinzirl et al. ( 2009)Q At high masses, 
all three models predict most galaxies are bulge-dominated. At low 
« 10 10 Mq) and intermediate (~ 10 10 - 10" Mq) masses, there 
is a dramatic difference: when the role of gas in suppressing disk 
destruction in mergers is accounted for, the model predicts the ex- 
istence of considerably more disk-dominated systems. 

The results are very similar when we adopt the semi-analytic 
approach from S08. We als o compare the semi-analytic model of 
Khochfar & Burkerj (2003), which relies on different cosmologi- 
cal assumptions and disk formation models, but still assumes the 
"threshold"-like (gas fraction-independent) model for bulge forma- 
tion in mergers, leading to a too bulge-dominated population. The 
authors attempt to address this with an age cut: enforcing that all 
mergers at z > 2 make no bulge. Figure|2]shows this is insufficient: 
many of these systems have too many mergers at late times to in- 
voke redshift-dependent dynamics or rapid subsequent cooling to 
"offset" the bulge formed. 

Figure [3] compares the median B/T (and scatter) as a 
function of stellar mass for each model to that observed. The 
predicted scatter owes mostly to differ ences in me r ger his - 
tory. We show the fits to observations i n [Balcells et al. (2007b), 
but others are similar jKochanek et alj [200 ll; lAllen et alj 120061 : 



iNoordermeer & van der Hulstl 120071 ; iDriver et alj 120071) . All the 
models predict some mass dependence, as merger histories are 
not completely mass-independent. However, at low masses, the 
gas fraction-independent models predict B/T ~ 0.4 — 0.5 (E/S0 
galaxies), much larger than observed. Note that both gas fraction- 
independent models give similar results; tuning the relative impor- 
tance of major/minor mergers does not resolve their problems. 

Figure [4] shows the results from the semi-analytic model for 
the same quantities considered in Figures □ &[2] The differences 

4 The observed B/T are determined in //-band; the models in Figuref2]pre- 
dict B/T in terms of stellar mass. In the semi-analytic model, the predicted 
star formation histories are used to predict the //-band B/T; generally we 
find this is slightly smaller than the stellar mass B/T, owing to the lower 
mass-to-light ratio in the younger stellar disk. 



between the predicted bulge-to-total distributions are dominated by 
the different treatment of the role of gas in mergers, rather than 
by the methodologies or different choices in the galaxy formation 
model. 

Figure [5] compares the predicted stellar mass functions in 
the semi-analytic model. Figure [6] shows the corresponding early- 
type fraction as a function of stellar mass, compared with the re- 
su lts from th e morphological ly classified stellar mass functions 
of iBell et alj j2003h (see a l so iKochanek et al.ll200ll; iBundv et al.l 
120051; IPannella et all 120061; iFranceschini et alj 12006. who obtain 
similar results). Although the gas-fraction independent models 
yield good agreement with the total stellar mass function, they pre- 
dict an excess of bulge-dominated galaxies at low masses. This is 
robust to the exact morphology or B/T cut used to define bulge or 
disk-dominated systems; here we adopt B/T ~ 0.4 as a dividing 
line a s this is a good p roxy, in simulations, for the concentration 
cutin lBelletai]d2003l) . 

We emphasize that we have incorporated the bulge formation 
prescriptions directly from H09 with the values of all related pa- 
rameters taken directly from the results of the hydrodynamic sim- 
ulations. We have not changed any of the other physical recipes in 
the SAM, nor have we re-tuned any of the baseline free parame- 
ters. We find it particularly convincing that the gas-fraction depen- 
dent bulge-formation prescription based on the simulations imme- 
diately produced excellent agreement with the observations at all 
mass scales, with no fine-tuning. 



3.2 Redshift Evolution 

Figure [7] repeats the comparison from FigureQ]at z = 1, 2, 3. Fig- 
ure [8] plots the median B/T and early-type fraction versus redshift 
and stellar mass. In the gas fraction-dependent model, higher gas 
fractions lead to more disk-dominated galaxies at high redshifts 
(see § 13. 3t . The evolution in g as fractions that we adopt (fitted to the 
observations in lErb et aL1l2006t) is relatively weak; if real evolution 
in gas fractions is stronger, it will further suppress B/T . 

Merger rates were higher in the past, so the gas fraction- 
independent models predict that galaxies should be more bulge- 
dominated. By z ~ 1 — 2, such models would predict that essentially 
all galaxies at all masses are bulge-dominated. 

At the highest redshifts, morphological determinations are dif- 
ficult; however, ther e are some m e asurem ents at z ~ 1 • We compare 
to the results f rom IBundv et alj j2005l) from the GOODS fields, 
IPannella et all fe006l) f rom t he GOODS South and FORS Deep 
fields, and lllbertetal] J2009I) from the COSMOS fields. Despite 
large uncertainties, the trend in the fraction of bulge-dominated 
galaxies with mass is qualitatively very similar to that observed 
at z — 0, but with a systematically lower bulge-dominated (E/S0) 
fraction at all masses. 

This agrees well with the expectations from the gas fraction- 
dependent model, and strongly conflicts with the gas fraction- 
independent model prediction. If anything, the data suggest that at 
high masses, bulges may be even more suppressed than the dissipa- 
tional model prediction; as noted above, this could be explained if 
gas fractions at these masses evolve more rapidly. Theoretical mod- 
els suggest, for example, that the efficient cooling (" cold mode") 
regime may persist to high m asses at high redshifts dKeres et"al] 
2005. 12009 : iDekel et alj|200^ . However, cosmic variance (not in- 
cluded in the quoted errors) remains a significant source of uncer- 
tainty at high masses. 

Figure|9]compares the redshift evolution of the predicted bulge 
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Figure 1. Predicted distribution of bulge-to-total stellar mass ratios in galaxies of different stellar masses at z = 0, from our semi-empirical model, with different 
models for bulge formation in mergers (see text). Left: The full (gas fraction-dependent) model from H09: bulge formation is suppressed in gas-rich mergers. 
Center: The "dissipationless" model: bulge formation is gas fraction-independent. Right: A simplified "threshold" model, also gas fraction-independent, with 
simplified treatment of the dependence on merger mass ratio. Accounting for the gas fraction-dependent physics of angular momentum loss in mergers, the 
predicted B/T is suppressed in low mass, gas-rich galaxies. 
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Figure 2. As Figure[TJ for galaxies with M t > 10 10 A/fry We c ompar e the 
models from Figure [TJ with observations from Weinzirl et al. 12009). We 
also compare the semi-analytic model from Khochfar & B urkerJ 42003b . 
which adopts different assumptions and does not allow mergers at z > 2 
to form bulges, but employs the simplified (threshold) model for bulge for- 
mation at z < 2. Only the gas fraction-dependent model matches the obser- 
vations; suppressing high-z bulges alone is insufficient to resolve the dis- 
crepancy in gas fraction-independent models. 



mass densit)|f] As expected from Figure [8] in the gas fraction- 
independent models, the predicted bulge mass density approaches 
the global stellar mass density at high-z. In gas fraction-dependent 
models, the trend is reversed: at high redshift, higher gas fractions 



5 Specifically, we plot the mass density in bulge-dominated galaxies, which 
is not the same as the absolute mass density in all bulges, but is closer to 
the observed quantity. At high redshifts z > 1.5 observed morphologies are 
ambiguous; we show the mass density in passively evolving red galaxies 
as a proxy. This may not be appropriate, but at £ < 1 the two correspond 
well, and the compactness, size, and kinema tics of the "passive" objects 
do appear distinct from star-forming objec t s iKrieket al. 2006 ; Toft et al. 
l2007MTruiillo et alj2007r , lFranx et al . 2008; Ge nzel et alj2008h . 
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Figure 3. Mean B/T as a function of galaxy stellar mass in observations and 
our semi-empirical models. Dot-dashed black lines show the ±lcr scatter 
in the full (gas fraction-dependent) model. The sha ded magenta range is the 
observed relation (with scatter) from Balcells et al. 1 2007b). The models are 
all similar at > V masses/luminosities, where disks are relatively gas-poor; 
at lower mass, the differences from Figures [Tl2l are apparent. 



suppress bulge formation. The latter trend agrees with the available 
observational results. 

Though important for bulges, we find that these differences 
have little effect on the global star formation rate density. "Quies- 
cent" (no n merger-driven) s tar formation dominates at all redshifts 
(see also IWolf et al.l I2005I: iMenanteau et al]|2006t iHopkins etal] 
l2006l : lNoeske et al.ll2007l ; |jogee et al.ll2009l) . Differences owing to 
e.g. systems "retaining" more gas for later consumption are much 
smaller than the uncertainties from e.g. star formation rate and stel- 
lar wind prescriptions. 



3.3 Why Does It Work? 

Why does the revised bulge formation model produce the observed 
trends? Figure [Tol illustrates the relevant physics. Consider the ob- 
served gas fractions of disks as a function of stellar mass, shown 
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Figure 4. As Figures [TJ & [2j but showing the results of the semi-analytic 
model, with either the full (gas fraction-dependent) bulge formation model 
or threshold (gas fraction-independent) model. Bottom: The results for 
M* > 10 10 Mq galaxies with B/T computed in the //-band, pre-selected 
to include only late-type galaxies (B/T < 0.75), compared to the observa- 
tions of Wei nzirl et al l 12009B . The semi-analytic and semi-empirical mod- 
els yield similar conclusions. 



in the Figure to be a strong decreasing function of mass|f| Now 



At z = 0, the ga s fractions sh own are based on measured atomic HI gas 
fractions; Bell & de Jong correct this to include both He and molecu- 

lar H 2 ; lMcGaughl l2005l) correct for He but not H 2 feannappanl 12004I) gives 
just the atomic HI gas fractions (this leads to slightly lower estimates, but 
still within the range of uncertainty plotted; H 2 may account for ~ 20 — 30% 
of the dynamical mass, per the measurements in ljogee et al.l2005h . We em- 
phasize that these gas fractions are lower limits (based on observed HI flux 



assume that every such galaxy doubles its mass in a merger or se- 
quence of mergerfl and consider the resulting B/T at each mass. 
This represents a major perturbation, and in any of the gas fraction- 
independent models would uniformly transform every galaxy into 
nearly 100% bulge, independent of mass. 

As discussed in § Q] however, gas loses angular momentum 
to internal torques from the stellar disk. If the disk is mostly gas, 
then the stellar bar/perturbation that does the majority of the work 
torquing the gas is relatively less massive and so less gas falls to 
the center and starbursts. In a collision of two pure gas disks, there 
is no stellar feature to torque the gas (in fact, accounting for typical 
orbital angular momenta, such a collision typically increases the 
disk specific angular momentum). 

Given these physics and the observed trend of gas fraction 
with mass, we expect low mass galaxies (M» < 10 9 Mq) to have 
just ~ 20% bulge even after a recent 1 : 1 merger. If mergers occur 
at higher redshift, when gas fractions are higher, the suppression 
is more efficient. It is possible for every galaxy to double its mass 
at z ~ 1 — 3 in major mergers and still lie below the local B/T- 
mass correlation (allowing some later mergers). In practice, many 
of these bulges are made by minor mergers, but the relative sup- 
pression is similar dWeinzirl et al.ll2009l : lHo"pkins et alj2009el) . 

This accounts for most of the observed trend in bulge frac- 
tion versus stellar mass. However, it is supplemented by two other 
trends (the origin of the weak scaling in even the gas fraction- 
independent models): more massive galaxies tend to have had more 
violent merger histories, and to have lower cooling rates (leading to 
less disk re-growth and lower gas fractions). 

FigureQTjjllustrates this in another way — we show the distri- 
bution of gas fractions in mergers as a function of merger redshift 
and stellar mass of the z = descendant galaxy (not instantaneous 
mass at the time of the merger). Gas fractions decrease with time 
overall, but low-mass galaxies are always more gas-rich than high- 
mass galaxies. This is for two reasons: (1) the Kennicutt law pre- 
dicts that star formation, and therefore gas consumption, is more 
efficient in high mass, higher surface density galaxies; (2) cooling 
and accretion are quenched in high-mass halos (in semi-analytic 
models, this is typically attributed to heating of "hot" quasi-static 
halos by AGN feedback). 

Figures [10] &QT] illustrate that the dividing line between gas- 
rich and gas-poor galaxies moves to higher masses at higher red- 
shifts, suppressing bulge formation even in massive galaxies at high 
redshifts. In semi-analytic models, this evolution occurs for three 
reasons: (1) cooling rates increase uniformly with redshift, (2) the 
critical halo mass dividing the "cold" and "hot" accretion modes 
(below which accretion is efficient) likewise increases with increas- 
ing cosmic gas densities at higher redshifts, and (3) AGN feedback 
becomes less efficient at high redshifts. For most of the z — popu- 



in some ann ulus). At z = 2, direct measurements are not available; the gas 
masses from lErb et all fe006l) are estimated indirectly based on the observed 
surface densities of star formation and assuming that t he z = Kennicut t 
law holds; other indirect estimates yield similar results JCresci et alj2009l) . 
7 This could be a single 1:1 merger, three 1:3 mergers, or even ten 1:10 
mergers; because the efficiency of bulge formation as calibrated in simu- 
lations is nearly lineal' in mass ratio, they give indistinguishable results for 
our purposes here, so long as they all occur within a short timescale (see e.g. 
i Naab & Burkertl2003l;lBournaud et al.l2005tlCox et alj2008l ; lHopkins et all 
2009b). In practice, this is unlikely at smaller mass ratios; cosmological cal- 
culations show that e.g. ten 1:10 merger s are a less likely channel for growth 
of spheroids than three 1:3 mergers jGottloberetal.1l200ll; iMaller et alj 
l2006tlZheng et alj 2007; Ste wart et al.l2008bl) . 
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Figure 5. Comparison of the z = mass function of all (left), bulge-dominated (center), and disk-dominated (right) galaxies. We compare the gas fraction- 
independent (thr eshold) and gas fraction-dependent (full) model, in the semi-analytic model, to the observed morphology-dependent mass functions from 
lBelletalj f2()03). Despite freedom in the SAM to adjust prescriptions, gas fraction-independent models overpredict the low-mass bulge-dominated population. 
Including the physical dependence on gas fraction reconciles this discrepancy. 
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Figure 6. Early-type (B/T > 0.4) fraction as a function of stellar mass. We 
compare with the z = o bserved trend from the morphologically classified 
mass functions in lBelletaljf2003h . 



lation, there is a relatively narrow "window" of efficient bulge for- 
mation, where cooling slows down, galaxies consume most of their 
gas, and mergers become efficient at forming bulges. 

Disk survival and bulge formation are at some level zero-sum. 
However, because various processes contribute to disk growth at 
different times, it is useful to ask whether disk "survival" or the 
suppression of bulge growth is more important. In order to answer 
this, we consider the following simple experiment: in each merger, 
we destroy as much of the disk gas and stellar mass as predicted by 
the dissipationless (gas fraction-independent) model; however, we 
only add to the bulge the mass that the dissipational (gas fraction- 
dependent) model says should be violently relaxed or lose angular 
momentum. The remaining mass is thrown away. This is not physi- 
cal, but it allows us to see whether our conclusions would change if 
bulge formation was suppressed per our full gas fraction-dependent 
model but disks were still relatively efficiently destroyed in mergers 
(per our gas fraction-independent model). 

We find that the results are almost unchanged from the full 
gas fraction-dependent model. The over-prediction in gas fraction- 
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Figure 7. As Figure[T] but evaluated at different redshifts z = 1 , z = 2, z = 3, 
for our gas fraction-dependent model (top) and gas fraction-independent 
model (bottom), implemented within the semi-empirical framework. The 
gas fraction-independent model predicts that the population is more bulge- 
dominated at higher redshifts (because merger rates are higher). Including 
the gas dependence reverses this: bulges are less dominant at high redshifts, 
as gas fractions are systematically higher. 

independent models of B/T is a problem of too much bulge for- 
mation, rather than too little disk survival from early times. This 
should not be surprising; various constraints suggest that low- 
mass galaxies are growing rapidly, and so disk masses are usu- 
ally dominated by recent accretion and star formation (see e.g. 
iBell & de Jon3l2000l : lGallazzi et al.ll20ol ; lNoeske et alj|200l and 
references therein). 



4 DISCUSSION 

A realistic model for how gas loses angular momentum and builds 
bulges in mergers qualitatively changes the expectations for the 
global distribution of bulge to disk ratios in the local Universe and 
the efficiency of bulge formation as a function of mass and red- 
shift. As demonstrated in simulations in H09, gas in mergers loses 
its angular momentum to internal torques from the induced non- 
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Figure 8. Early-type fraction at different redshift s, in the semi-empi rical 
mode l. Green points sh ow the obser vations from iBundv et alj J2005l cir- 
cles), HjbertTeTal] (2003, squares), and Pannella et al. 12006, triangles) from 
the COSMOS, DEEP and GOODS fields at z « 1 (compare the green dot- 
dashed line). Ignoring the gas-fraction dependence predicts that even low- 
mass populations are > 80% elliptical at z > 1. 
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Figure 9. Integrated mass density in bulge-dominated galaxies in the 
semi-empirical models as a function of redshift, compared to ob- 
servations (points). Contr ast the observed total stellar mass density 
I Hopkins & Beacoml l2006t dot-dashed line). Observations are compiled 
from the morpholo gically-selected sam ples of iBell et al. fe003l black X ), 
IBundv et alj_ 20051 12006L red circles), I Abraham et afl~ 200l violet dia- 
monds), andlDaddi et al J fe005l c yan square), and color-sele cted samples of 
Labbe et al. (2005, green square). Ivan 

Dokkum et allfeOOd orange square), 
and lGrazian et alJ i2007. magenta stars). 
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Figure 10. Top Left: Obser ved median gas fracti on as a func tion of stel- 
lar ma ss in z = dis ks. fromlBell & de Jond 1200 ll diamonds), [Kannappan 
120041 squares), and McGaugh 120051 circles), with an analytic fit (solid 
line) and ±lcr scatter (dashed). Dotted (blue) line compares median gas 
fractions in the semi-analytic model. Top Right: Predicted bulge-to-total ra- 
tio in the remnants if every such galaxy immediately underwent a series of 
mergers that doubled its stellar mass (style as Figure|5J. Lines correspond to 
the different models: full/gas fraction-dependent (solid black), dissipation- 
less/gas fraction-independent (dashed blue), and simplified/threshold (dot- 
ted red). Compare the observed trend (shaded). Bottom: Same, but based on 
the gas fractions at z ~ 2 lErb et al. 2006). Accounting for the full model 
results from simulations, high gas fractions in low-mass systems suppress 
bulge formation. 
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Figure 11. Galaxy gas fractions as a function of the z = stellar mass of the 
descendant and the time at which the merger occurred (in the semi-analytic 
model). Specifically, considering galaxies which have the given z. = stellar 
mass, and whose main branch progenitor experienced a merger at the given 
redshift (scaled here in equal-cosmic time units), we plot the median gas 
fraction of the progenitors just before the merger(s). All galaxies were more 
gas-rich at high redshift, but low-mass galaxies are more gas-rich than high 
mass galaxies (on average) at all epochs. 
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axisymmetry in the stellar distribution in the same disk as the gas. 
Therefore, in the limit of very gas rich mergers (with little or no 
stellar mass in the disks to do such torquing), angular momentum 
loss and destruction of the gas disk is inefficient (to lowest order, 
the efficiency scales oc (1 — / gas )). 

We incorporate this idea into a cosmological framework us- 
ing two different approaches: a simple semi-empirical model and 
a full semi-analytic model. The conclusions appear robust to this 
choice; we have also experimented with a number of variations to 
the semi-empirical model (detailed comparison of the resulting dif- 
ferences for merger rates will be the subject of future work; we 
briefly summarize the variations in Appen dix IO. and find similar 
results in each case. IStewart et al.l J2009h also consider an inde- 
pendent semi-empirical model, using different halo merger trees, 
observational constraints, and more qualitative prescriptions for 
behavior in mergers (rather than the detailed simulation results 
here), and reach the same conclusions. It appears that, regardless 
of the methodology, this picture of morphological transformation 
in mergers makes qualitatively different predictions from the "tra- 
ditional" picture which ignores the dependence on gas fraction. 
Because, in the local Universe, lower mass disks have large gas 
fractions (/ gas > 0.5 at M* < 10'° Mq), this in particular implies 
stronger suppression of bulge formation at low masses. Even if the 
dissipationless (halo) merger histories were self-similar across dif- 
ferent mass scales, this trend of gas fraction with mass breaks this 
self-similarity and largely explains the observed trend of bulge frac- 
tion with mass. 

Other factors are expected to weakly break this self-similarity 
(the non-linear dependence of galaxy baryonic mass on halo mass, 
and weak scale dependence in the dark matter merger histories), 
but simulations and empirical constraints have found that account- 
ing for these factors alone predicts too weak a trend. On the other 
hand, allowing for the role of gas-richness (implicitly related to the 
nature of cooling and star formation efficiency) provides a natu- 
ral and predictive explanation for the observed mass-morphology 
relation. 

Moreover, whereas simulations without large gas supplies, or 
gas fraction-independent models have difficulty reproducing the 
abundance of "bulge-less" or low B/T disks (B/T < 0.1 —0.2), 
this problem is largely solved by incorporating the gas fraction- 
dependent bulge formation prescription. In practice, if galaxies 
have sufficiently low star formation efficiencies in order to retain 
very large gas fractions, as implied by observations of these sys- 
tems, and especially if the gas distribution is very extended (as in 
e.g. low surface brightness systems), then it is possible to produce 
systems with such low B/T even in the immediate aftermath of 
a major merger, let alone in systems which have gone ~ lOGyr 
since the last such merger with significant subsequent accretion. In 
practice, many bulges of low (and, to a lesser extent, intermediate) 
mass systems are predominantly made by minor mergers, making 
the required gas-richness even more moderate. 

At low redshifts, this is primarily relevant for small-mass 
galaxies, where gas fractions are high. However, observations in- 
dicate that gas fractions should be larger at high redshifts, such 
that even some ~ L* galaxies appear to have / gas > 0.5 at redshifts 
z > 2 — 3. At increasing redshift, therefore, the effects described 
above propagate to higher and higher masses, and become impor- 
tant for an increasingly large proportion of the galaxy population. 

This is seen as well i n recent cosmological simulations 
dGovernato etai]|2007l. l200i) : the presence of gas leads to non- 
destructive mergers and rapid disk re-assembly. The role of stellar 
feedback is to preserve gas, preventing it from all turning into stars 



at early times or in very small galaxies, and to remove low an- 
gular momentum gas from the disk (keeping the gas at relatively 
large radii, comparable to observed disk sizes). Even at z ~ 0.5, 
disks experience major mergers in these simulations and rapidly re- 
build. These processes, it appears, are now being seen in a growing 
number of observations dHammer et alj|2005l: IZheng et alj 12005 ; 
i Truiillo & PohlerJ 120051: iFlores et al.l 120061 ; IPuech et alj l2007li . 
2008; Atki nson et al.1 12007). The suppression of bulge formation 
in the merger owing to its gas-richness, combined with rapid large- 
scale gas inflows concurrent with and subsequent to the merger, 
both contribute to the rapid re-building of a disk-dominated galaxy 
follow ing the event ( Springel & Hernquist 2005,; iRobertson et al.l 
l2006al) . The inclusion of this continuous accretion means that a 
very large gas fraction / gas > 0.5 is not needed to build a disk- 
dominated galaxy (as would be the case if the system received no 
new gas after the merger); the more moderate bulge suppression 
that would result from intermediate gas fractions is sufficient. 

At higher redshifts, merger rates increase, but in gas fraction- 
dependent models this is more than offset by the increase in disk 
gas content and cooling efficiencies, and thus low B/T ratios are 
maintained. However, because the rates of these processes (merg- 
ers and cooling/accretion) scale with the Hubble time while the in- 
ternal relaxation times of disks are relatively constant, disks should 
be more turbulent and the relaxation from the building/merger/re- 
formation process should be evident in more systems, as thicker, 
more dispersion-dominated disks with enhanced incide nce of struc- 
ture and dumpiness. In fact, this is widely observed, l Weiner et al.1 



l2006l;lAtkinson et alj2007; Puec h et alj2007bl : 



Yang et al. 2008) especially approaching z ~ 2 (Bouche et al. 200' 



Neichel et alBoQ: 



Shapi ro et alj|2008l; Ivan Starkenburg et al.ll2008l) , and agrees well 
with the structure of simulated disks underg oing rapid simultane- 
ous accretion and post-merg er reformation jRobertson & Bullockl 
l2008bl ; lGovernato et alj|2008l) . 

Other observational signatures of this process include "in situ" 
bulge formation: rather than the conventional wisdom, in which 
most classical (merger-induced) bulges are formed by mergers that 
destroyed the progenitor disk at some early time, with a later re- 
accreted disk around them, the dissipational models predict that it 
is much more likely, especially in gas-rich (low-mass) systems, that 
such bulges were formed by more recent mergers (many of them 
minor) that formed a small bulge out of a large disk, without de- 
stroying most of the disk. As such, the properties of the stellar pop- 
ulations in the bulge should be reasonably continuous with those in 
the disk. Note that, in an integrated sense, the bulge may have older 
or more metal-rich stars, but this is true even in a pure disk system 
if one compares the central regions to the outer regions; moreover, 
star formation may continue in the disk after the merger, whereas it 
may not in the bulge. But, in terms of e.g. the full stellar population 
gradients and even resolved color gradients, there should be conti- 
nuity between disk and bulge, with properties such as their metal- 
licities correlated. Indeed, high resolution resolved spectroscopic 
observations of disk+bulge systems have begun to show that such 
continuity is the norm, i n metallicity and q-enrichment, stellar pop- 
ulation ages, and colors i Moorthv & Holtzmarl2006|; Peletier et alj 
20071; iDomfnguez-Palmero & Balcellsl 120081; rMorelli et alj|2008l 
MacArthur et alj2008l) .~ 

Observations and simulations have shown that the two- 
component (gas+stellar) nature of mergers leads to a two- 
component structure of bulges, with a low-density outer compo- 
nent reflecting disk stars that are violently relaxed in the merger and 
dens e inner component built b y the gas that loses ang ular momen- 
tum jMihos & Hernqui st! 1994b . In a series of papers, Hopki mTet al.1 
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( 2008bl. l20093ldLl2008al . [2009ch demonstrated that with sufficiently 
high-dynamic range observations of spheroid structure and/or stel- 
lar populations, these can be decomposed, yielding empirical con- 
straints on the efficiency of angular momentum loss in gas. Con- 
sideri ng a large sample of ellipticals and recent merger rem- 
nants teothberg & JosephfeOoHlLauer et alj2007l ; lKormendv et al.l 
2009), the mass fraction of spheroids in the dissipational com- 
ponent appears, to lowest order, to tra ce typical disk gas fra c- 
tions as a function of mass; however, iHopkins et al.l 1 2009al lch 
note that at low masses, this appears to asymptote to a maxi- 
mum dissipational/starburst fraction ~ 0.3 — 0.4. Independent tests 
dBalcells et alj2007 d ; lHopkins et alj2008bl;lKormendv et al.l2009T) 
and constraints from stellar populations {McDermid et al. 2006; 
ISanchez-Blazquez et alJl2007l ; lReda et alj|2007h yield similar con- 
clusions. If angular momentum loss in mergers were always effi- 
cient, the natural expectation would be that this should rise at low 
galaxy masses as the progenitor disk gas fractions, i.e. ~ / gas — > 1. 
If, however, angular momentum loss becomes less efficient in these 
extremely gas-rich systems, the observed trend is the natural pre- 
diction: if the fraction of gas losing angular momentum scales as 
adopted here, then the dissipational fraction of the bulge formed 
from disks with gas fraction / gas is not ~ / gas , but ~ /gas/(l +/gas)> 
i.e. asymptoting to the values observed for all / gas ~ 0.5 — 0.9. Ex- 
panded samples with high-resolution photometry, kinematics, and 
stellar population constraints would provide powerful constraints 
to extend these comparisons. 
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APPENDIX A: DETAILS OF THE HOPKINS ET AL. 2009 
(H09) MODEL 

The full gas fraction-dependent ("new") model for the efficiency 
of bulge formation as a function of gas fraction and merger prop- 
erties, discussed in the text, is derived in detail and compared to a 
large suite of hydrodynamic merger simulations spanning the rel- 
evant parameter space in H09. We refer the reader to that paper 
for a detailed discussion of the physical motivation for the model 
and numerical tests. Here, we briefly present the key physically- 
motivated prescriptions which are shown therein to provide a good 
fit to the simulation behavior and provide our implementation of 
the hydrodynamic results in analytic or semi-analytic cosmological 
models. 

Define: m»,di s k the mass of stars in the disk component and 
m*, bulge the mass of stars in the bulge component, with m, — 
"Vdisk + bulge the total stellar mass; m co id the mass of cold gas 



in the galaxy (assumed to all reside in the disk). Here, mbar = 
m» + m co id is the baryonic mass of the galaxy, and r<i is the effective 
radius of the disk. The disk gas fraction is then / gas = m co id / (m co id + 
m*,disk) and the disk mass fraction is ./disk = (mcoid + m*.disk)/mbar[3 
The model outlined in § 4.3 of H09 then gives a simple set of 
formulae to compute the fraction of cold gas in both merger progen- 
itors that will participate in a central starburst (as a consequence of 
losing its angular momentum), and as a result be removed from the 
disk and added to the bulge component of the remnant. The mass 
of gas that participates in the starburst is 



«cold (< fcrit), 



(Al) 



or the mass of cold gas within a critical radius r c ,it where angular 
momentum loss from the gas disk in the merger is efficient. 

H09 derive the following relation for r tT H (their Equation (7)): 



'"oil 



= a{\-f m )fa*F{0,b)G{p) 



where 



i 



3/2 



1 



+ [b/r d ] 2 J i-n D /n d 

contains the dependence on the merger orbit, and 

2 1 1 



G( M ) 



i+h 



(A2) 



(A3) 



(A4) 



contains the dependence on the merger mass ratio n = nvz/mi. 
There is some ambiguity in the most relevant definition of "mass 
ratio" (e.g. lStewartll2009l) : because we are interested in the gravita- 
tional dynamics induced by passage/merger, we find the best match 
to behavior in simulations by defining the masses of interest (mi 
and mo) as the baryonic plus tightly bound (central) dark matter 
(specifically, dark matter mass within a core radius of r s , the NFW 
scale length = r yn /c, roughly ~ 10% of r v h or ~a few disk scale 
lengths; the results are not strongly sensitive to the choice within a 
factor of ~ 2 — 3 in r). This is important, in particular for low-mass 
galaxies which are observed to be dark-matter dominated (and have 
high /g ;l s); simply considering the stellar or baryonic mass ratio sig- 
nificantly under-estimates the "damage" done to the disk relative to 
the results of high-resolution simulations. 

In Equation I A3 1 above, b is the distance of pericentric pas- 
sage on the relevant final passage before coalescence (< a couple 
rd for typical cosmological mergers), but H09 show that for sys- 
tems which will actually merge, as opposed to fly-by encounters, 
identical results are obtained in the limit b — > 0, since the final coa- 
lescence — regardless of initial impact parameter — will be largely 
radial. The quantity fi is the orbital frequency at pericentric pas- 
sage and Qd is that of the disk; 



^'peri » d 



cos (0) 



x/2(1+ai)[1 + (b/ ri ) 2 r 3/4 cos (9), 



(A5) 



where 9 is the inclination of the orbit relative to the disk (v pe ,i is the 
velocity at pericentric passage, v c the circular velocity at the im- 
pact parameter b). We draw 9 randomly for each merger assuming 



8 In simulations, the following descriptions are most robust if these quanti- 
ties are denned just before the final coalescence/merger where the gas will 
lose angular momentum and stars will violently relax. Defining them at ear- 
lier times means that other quantities such as the star formation efficiency 
will be important for determining how much gas is actually present at the 
time of interest (the merger). 
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disks are isotropically oriented (although this affects primarily the 
scatter, not the mean relations predicted). 

Adopting b = 2rd, typical of cosmological mergers, or taking 
the limit b — > for final radial infall gives a similar scaling: 

aF(9,b) = — . (A6) 

v ; 1- 0.42 ,/TT7I cos (9) 

The parameter a in Equation l lA2t subsumes details of the stellar 
profile shape and bar/driven distortion dynamics in a merger. The 
precise value above is obtained from the numerical simulation re- 
sults in H09, but is also derived analytically therein. 

The disks are assumed to follow an exponential profile, so the 
fraction of the cold gas mass within rent is given by 

/burst = = 1 - (1 + r cn ,/r d ) exp (~r cnt /r<i). (A7) 

Generally, evaluating these quantities with the cold gas and stellar 
content being the sum of the two merging galaxies is a reasonable 
approximation to the numerical results (representing the fact that 
the burst takes place near or shortly after the time of coalescence); it 
is also straightforward to independently apply the equations above 
to both primary and secondary at the time of the merger. In practice, 
it makes little difference (even assuming that the entire gas mass of 
the secondary is part of the burst is a good approximation for all 
but the most rare very major mergers with mass ratios close to 1:1). 

Following H09 § 4.3.3, we further assume that the mass of the 
stellar disk in the primary (more massive) galaxy that is transferred 
to the bulge/spheroidal component from the disk and violently re- 
laxed is 

»Idisk, destroyed = /J m*. disk- (A8) 

Furthermore, the whole stellar mass of the secondary, m%, is also 
added to the spheroidal remnant. A more detailed accounting of e.g. 
the efficiency of viole nt relaxation as a function of radius (see H09, 
Hopkins et al. 2008c) in the stellar disk suggests that one should 
further multiply the right-hand side of Equation lA8l bv a factor ( 1 + 
fl 3 ) _1 ; in practice including this weak dependence on fi makes 
no difference. 

As described in the text, the gas fraction-independent or "dis- 
sipationless" model treats gas in the same way as stars in terms of 
the efficiency of merger-induced bulge formation. Gas which is as- 
sumed to be transferred to the spheroid in the merger is still allowed 
to lose angular momentum and participate in a central starburst as 
opposed to being violently relaxed (so the system is not, formally 
speaking, dissipationless), but the gas mass fraction that does so 
scales with merger properties in the same manner as the dissipa- 
tionless violently relaxed stellar fraction. This amounts to replacing 
our more physical model for the gas behavior with the analogue of 
Equation |A8l the gas mass which goes into a starburst in the merger 
is simply OTburst = M»i c( ,ki, independent of gas fraction. Since the be- 
havior of interest is the dependence of bulge formation efficiency 
on gas fraction, we can instead still choose to adopt Equation lA2l for 
the behavior of the gas, but omit the dependence on disk gas frac- 
tion (the (1 — /gas) term), equivalent to assuming that the gas always 
loses angular momentum as efficiently as in a relatively gas-poor 
merger (formally computing Equation lA2l with / gas — > 0, regardless 
of the actual disk gas fraction). This choice makes no difference 
in our calculations and comparisons throughout the paper: the de- 
pendence on e.g. orbital parameters in Equation lA2l onlv affects the 
scatter in B/T, not the mean trends, and the scaling with mass ra- 
tio is very similar. The important difference between the models is 
entirely contained in the (1 — / gas ) scaling. 

We emphasize that there are no free or "tuned" parameters 



in this model. All of the parameters and values given above are 
determined in H09 based on physical models and hydrodynamic 
merger simulations, and are adopted without adjustment in both 
the empirical and the semi-analytic model. 



APPENDIX B: DETAILS OF THE S08 SEMI- ANALYTIC 
MODEL 

Here, we give further details of the semi-analytic model in which 
we embed the above models for bulge formation. Thi s mode l 
has been described in detail in ISomerville & Primackl d 1999b : 
ISomerville etal] fcOOll) a nd a version with major u pdates has re- 
cently been described in ISomerville eTaT] j2008aL S08). We re- 
fer the reader to those papers for details, but briefly summarize 
the main ingredients of the models here. The models are based 
on Monte Carlo realizations of dark matter halo merger histories, 
constructed using a slightly modified version of the method of 
ISomerville & Kolattl Jl999h . as described in S08. The merging of 
sub-halos within virialized dark matter halos is modeled by com- 
puting the time required for the subhalo to lose its angular momen- 
tum via dynamical friction and fall to the center of the parent halo. 
The S08 model takes into account the mass loss of the orbiting dark 
matter halos owing to tidal stripping, as well as tidal destruction of 
satellite galaxies. 

The model discriminates between cooling and accretion in two 
different regimes: "cold mode", in which the gas cooling time is 
less than the free-fall time, and "hot mode", in which the opposite 
condition holds. In the cold mode, gas accretion onto the disk is 
limited only by the rate at which gas is incorporated into the virial- 
ized halo. In the hot mode, the rate at which gas can lose its energy 
and cool via atomic cooling is estimated. This cooling gas is as- 
sumed to initially settle into a thin exponential disk, supported by 
its angular momentum. Given the halo's concentration parameter, 
spin parameter, and the fraction of baryons in the disc, angular mo- 
mentum conservation arguments are use d to compute the scale ra - 
dius of the exponential disc after collapse ( ISomerville et alj2 008b). 
The models make use of fairly standard recipes for the suppression 
of gas infall owing to a photo-ionizing background, star formation, 
and supernova feedback, as described in S08. Quiescent star for- 
mation, in isolated dis cs, is mod eled using the empirical Kenni- 
cutt star-formation law l lKennicutll 19981) . A critical surface density 
threshold is adopted, and only gas lying at surface densities above 
this value is available for star formation. Cold gas may be reheated 
and ejected from the galaxy, and possibly from the halo, by super- 
nova feedback. Heavy elements are produced with each generation 
of stars, and enrich both the cold and hot gas. 

Galaxy mergers trigger the accretion of gas onto supermassive 
black holes in galactic nuclei. Energy radiated by black holes dur- 
ing this 'bright', quasar-like mode can drive galactic-scale winds, 
clearing cold gas from the post-merger remnants. In addition to 
the rapid growth of BHs in the merger-fueled, radiatively efficient 
'bright mode', BHs also experience a low-Eddington-ratio, radia- 
tively inefficient mode of growth associated with efficient produc- 
tion of radio jets that can heat gas in a quasi-hydrostatic hot halo. 
Thus, once the black hole has grown sufficiently massive, this "ra- 
dio mode" feedback can completely switch off further cooling, 
eventually quenching star formation in massive galaxies. Note that 
the radio mode feedback is assumed only to be effective in heating 
gas that is being accreted in the "hot mode" described above. The 
semi-analytic models contain a number of free parameters, which 
are summarized in Table 1 of S08. These parameters were chosen 
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in order to reproduce certain observations of nearby galaxies, as 
described in detail in S08 (most importantly, the stellar mass func- 
tion, cold gas fraction of spirals as a function of stellar mass, and 
the metallicity-mass relation zero-point). In this work we use the 
fiducial values of all parameters, as given in S08. 



galaxy gas fractions. Those are the parameters that effect the model 
conclusions; to the extent that the variations above provide a means 
to interpolate between the empirical constraints on these quantities, 
they yield similar conclusions. 



APPENDIX C: VARIATIONS OF THE SEMI-EMPIRICAL 
MODEL 

An advantage of the semi-empirical model approach is that it is 
trivial to alter various assumptions, provided that the key observa- 
tional constraints (e.g. galaxy masses and gas fractions, in this case) 
remain satisfied. We have experimented with a number of such vari- 
ations, and find that the conclusions in this paper are robust. A 
detailed discussion of these variations and their consequences for 
merge r rates and galaxy me rger histories will be presented in future 
work dHopkins et al.l20 09e). Here, we briefly outline some of these 
choices. 

Changing the adopted dark matter merger rates makes lit- 
tle difference. We have re-calculated ou r model adopting the halo 
merger trees from lStewart et aL (2008a), using the subhalo-based 



merger trees from IWetzel et alj (2003), ar, d implementing vari- 
ous approximations for a "merger timescale" (delay between halo- 
halo and galaxy-galaxy merger, given e.g. by some approximation 
to t he dynamical fricti on time) from Bovlan-Kolchin et al. (2008) 
and I Jiang et ail f 2008). We have also re-calculated the model us- 
ing different observational constraints to determine the galaxy 
masses populated in halos or su bhalos of a given m ass; using 
the monotonic-ranking method of IConrov et alj J2006|) applied to 



the re dshi ft-dependent galaxy mass fu nctions from Fontana et alj 
(2006) or Perez -Gonzalez et al. (2008), or using the z = con- 
straints determine d directly from clustering of SDSS galaxies in 
IWang et al.1 120061) and ignoring redshift evolution in this quantity. 
We find in each case a similar result. 

We can also change the methodology, to essentially a toy 
model (between our strict semi-empirical approach and a semi- 
analytic model) fitted to observations. We assume at each timestep 
that all galaxies in halos below 10 12 Mq simply accrete cold gas 
proportional to their halo growth (this corresponding to what is 
seen in simulations and observations; low mass galaxies accrete 
efficiently in the "cold mode" with gas coming directly into disks 
on a dynamical time, h igh mass galaxies are "quenched"; see e.g. 
iKeres etalj|2005t l2009h . We can then determine a star forma- 
tion rate either by forcing the galaxies to move along the em- 
pirically constrained M ga i(Mh a io) track (our default method, based 
on the observational constraints described above), or by assum- 
ing the galaxies lie on a fixed size-mass relation and using the 
iKennicuttl fl998l) relation (both appear to evolve onl y weakly with 



redshift; see 



Bouc heetal.il2007l; IPuech et al.1 120081 ; iGenzel et al.1 



l2008l ; lvan Starkenburg et 311200% or by assuming th ey live on the 
fitted SFR-galaxy mass relation in lNoeske et alj ( l2007l) . We can fur- 
ther allow some "outflow," with an efficiency M out a ow oc M, fitted so 
as to yield a good match to the /gas-galaxy mass relation observed 



as to yield a good match to the / aas -galaxy mass relation observed 
at different redshifts (e.g. iBell & de Jondl200d;lMcGaughl 120051 ; 



IShaplev et al]|2005l ; lErb et alJl2006l ; ICalura et al.l(2008l) . Recall, in 
our "default" semi-empirical model, we simply assume whatever 
inflow/outflow necessary such that the galaxies always live on the 
/g as -galaxy mass relation fitted to these observations; but appropri- 
ately tuned, this simple toy model can yield a similar result. 

Ultimately, we find that these choices make little difference, 
so long as they yield the same distribution of galaxy masses and 
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